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Abstract. This paper is devoted to the muhigrid convergence analysis for the hnear sys- 
tems arising from the conforming hnear finite element discretization of the second order 
elliptic equations with anisotropic diffusion. The multigrid convergence behavior is known 
to strongly depend on whether the discretization grid is aligned or non-aligned with the 
anisotropic direction and analyses in the paper will be mainly focused on two-level algo- 
rithms. For an aligned grid case, a lower bound is given for point-wise smoother which 
shows deterioration of convergence rate. In both aligned and non-aligned cases we show 
that for a specially designed block smoother the convergence is uniform with respect to 
both anisotropy ratio and mesh size in the energy norm. The analysis is complemented with 
numerical experiments which confirm the theoretical results. 



1. Introduction 

In this paper we will study multilevel methods for anisotropic partial differential equations 
(PDEs) by finite element (FE) methods and in particular we will analyze the convergence 
behavior of some methods for anisotropic diffusion equations on grids that are either aligned 
or non-aligned with the anisotropy. 

There are already many convergence results in the literature for multilevel method on 
anisotropic problems when the underlying FE grid is aligned with the anisotropy direc- 
tion. The case of constant anisotropy was considered by Stevenson [D, [2], who established 
a uniform convergence of the V-cycle multigrid methods. Main tools in his analysis are 
the so called classical smoothing and approximation properties (see Hackbusch [3]). The 
case of "mildly" varying anisotropy was analyzed in a work by Bramble and Zhang [1] . Us- 
ing a different theoretical framework developed in Bramble, Pasciak, Wang and Xu [5] and 
Xu |6], Neuss [7] also showed uniform convergence of the V-cycle algorithm for anisotropic 
diffusion problem. More recently, Wu, Chen, Xie and Xu |8j analyzed V-cycle multigrid 
with line smoother and standard coarsening, and V-cycle multigrid with point Gauss-Seidel 
smoother and semi-coarsening at the same time, and they were able to prove convergence 
under weaker assumptions on the regularity of the solution of the underlying PDE. An- 
other technique, based on tensor product type subspace splittings and a semi-coarsening was 
proposed and analyzed by Griebel and Oswald [9j. They have shown uniform and optimal 
condition number bounds for multilevel additive preconditioners. 
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The aforementioned theoretical convergence results on multigrid methods for the anisotropic 
diffusion equations are, however, all carried out under one main assumption that the anisotropy 
direction is aligned with the mesh. But such an assumption is not always satisfied in practice. 
In this paper, we will make an attempt to develop uniform convergence theory in certain 
cases when this aligned grid assumption is not satisfied. More specifically, we will study the 



problem (2.1) in the case that f2 is a square domain triangulated by uniform grids that is 



rotated by an angle uj G [0,7r]. We have grids that are not aligned with anisotropy except 
for special cases that Ci; = 0, Ci; = |ora; = ^. 

For this special class of domains and grids, we will design a two-level method and prove 
its uniform convergence (with respect to both anisotropy and mesh size). We are not yet 
able to extend our theoretical analysis neither to multilevel (more than two levels) case, nor 
to more general anisotropic problems. We hope however that the analysis presented here, 
even though in a special case, can be extended to handle more general anisotropic problems. 

We would like to point out that our work was partially motivated by some recent theoretical 
results for nearly singular problems (see p^). Indeed, anisotropic diffusion equation gets 
more nearly singular when anisotropy gets smaller. Techniques such as line smoother or semi- 
coarsening correspond in some way to space splittings of the so called near kernel components 
of the anisotropic diffusion problem. We refer to pU] for description of such splittings. 

The rest of this paper is organized as follows. In Section |2] we introduce the notation 
and preliminaries. In Section [3] we state the main result. We then prove several stability 
and interpolation estimates for coarse grid interpolant in Section |4] and for the fine grid 
interpolant in Section [5j In section [6] we prove the main theorem, already stated in section [3j 
Numerical experiments, which verify the theory are given in section [7] 

2. Preliminaries and notation 
Consider the anisotropic diffusion equation on a square domain i7 C M^: 

'^xx ^^yy fi ^1 



^^■^^ ^ M = 0, on dQ, 

where e > is a constant. We are interested in the case when e — )■ 0. The weak formulation 



of (|2l|) is: Find u e H^{fl) such that 

(2.2) a{u,v) = {f,v), yveH',{n), 

where 



a{u, v) = {dxudxV + edyudyv)dxdy, and (/, v) = {fv)dxdy. 
Jn Jq 
We consider family of computational domains Q obtained by rotations of a fixed domain 
Qq = (—1, 1)^ around the origin. The angles of rotation are denoted by u and we consider u G 
[0,7r], since this covers all the possible cases of alignment (non-alignment) of the anisotropy 
and the FE grid. 

We assume that we have initial triangulation To of the domain Qq, obtained by dividing 
f^o into N X N equal squares and then dividing every square into two triangles. Then we 
rotate Qq around the origin to obtain the computational domain Q and its triangulation Th- 
The finite element function space associated with Q and Th will be the space of piece-wise 
continuous linear functions with respect to Th and we denote this space by Vh- One may see 



three such domains shown in Figure 2.1 In such setting one case of grid aligned anisotropy 



corresponds to f2 = Hq, or equivalently, u = 0. 



TWO-LEVEL METHOD FOR ANISOTROPIC DIFFUSION EQUATIONS 



3 




Given a coarse mesh Th-, assume the fine mesh Th is obtained from Th by sphtting each 
of the triangles in the triangulation Tff into four congruent triangles. One clearly then 
has h = H/2. The spaces of continuous piece- wise linear functions corresponding to the 
partitions Th and Th are denoted by Vh and Vh- As it is customary 1^, Ir will denote the 
nodal interpolation operators mapping to I4 and Vh respectively. 

For the analysis of the two-level method, we introduce partition of unity where 
L — r ^'"°'^~^'"'" 1 , with Umax = max {y}, y^nin = mill {y}- We define 9i as follows: 

r {t-l)h<y-y^in<th, 

(2-3) Oi{y)=< zh<y-ymin<{t + l)h, 

[ 0, other. 

Note that each 6i is piece-wise linear in the y variable and is constant in x. Moreover, each 
9i is supported in the i-th strip (i — l)h < y — ymin < + 1)^- 

Denote the set of triangles in Th including nodes in the i-th strip by 71, and Let 

= U r, 

reTi 

Vi = {v e Vh, supp V C Qi}, 

then the two-level method with line smoother Vi{l < i < L) and coarse grid Vh can be 
written as 

L 
i=l 

Let be a triangle with vertices {(xj, which we assume ordered counter-clockwise. 

For a given edge E e dK, with 

E = {{xi,yi), {xj,yj)), j = 1 + mod(i,3), 

we denote 

(2.4) S§y^^^{y,-y,), 

we also denote with {xE,yE) the coordinates of the vertex of K which is opposite to E. In 
another word, if ii^ = {{xi,yi), {xj,yj)) then {xE,yE) = {xk,yk), where k i and k 7^ j. Let 
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w be a linear function on K. If we set = v{xE,yE), then, it is easy to check that 



)4- 



K 



EedK 



3. Convergence of the two-level method 

We will first prove that even in the aligned case the most common point-wise smoothers 
will result in two-level method whose convergence deteriorates when e tends to zero in equa- 



tion (2.1). The result is as follows. 



Theorem 3.1 (Lower bound and deterioration of the convergence rate). In case of grid 
aligned anisotropy (i.e. u = 0, the energy norm of the error propagation operator 

corresponding to the two-level iteration with coarse space Vh and point-wise Gauss-Seidel 
smoother can he hounded helow as follows: 

(3.1) \\ETL\\l>l-C{e + h^), 

with constant C independent of e and h. 

This result follows from the following two-level convergence identity (proof can be found 
in pTl Lemma 2.3]): 

Lemma 3.2. The following relation holds for the two-level error propagation operator Etl = 
{I~T){I-Ph): 

— sup 



(3.2) ll-S^Tilla = 1 - where if = sup .. 

K „pv WvWi 



where \\v\\1 = inf ^ ll^^fclla ^'^^ ^* ^■^ ■)*- orthogonal projection on Vh. 

T.t-"i=^ k 



Proof of Theorem 3.1. From the Lemma 3.2 we can immediately see that to prove the esti- 



mate (3.1) we need to show that 



(3.3) A- up 11'^- "•>""•> 



here the quantity K is the same as in Lemma 3.2 From the proof of pTl Theorem 4.5], we 
also know that 

(3.4) K > sup 

where Qh is the (■, ■) orthogonal projection on Vh- In the case of angle of rotation w = 0, the 
computational domain is f2 = = (~1) 1)^- We assume h = 1/n (n is even) and consider 
the 2n x 2n partition with vertices {xj,yk), Xj = jh and y^ = kh, j = — n, ■ ■ ■ ,n,k = 
—n, ■ ■ ■ ,n, then the corresponding coarse gird is with vertices {x2j,y2k), j = ~n/2, ■ ■ ■ ,n/2, 
k = -n/2,-- - ,n/2. 

For any given v & Vh, since 



-n k=—n 
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and 

n/2 n/2 
j=-n/2 k=-n/2 

then the interpolation Ih is stable in the L2 norm, i.e. H/^'i'll ^ H'^^H, for all v G Vh, provided 
that ^ < 1. Now consider a function t>o G Vh, supported in the closure of (—1, 1) x (0, 2/i) 
and defined as 

Vo{xj,yi) = Vo{xj,h) = 1 - \j\hj = -n, ■ ■ ■ ,n, 
and fo is at any other vertex. Note that 

IhVq = 0. 

From the stability of Ih in the L2 norm, which we have just shown we get: 

\\lHiI-QH)vof < UI-QhHW, 

\\(i-i^)(i-Q^)vor < ui-QnHr. 

Using these estimates and the fact that IhQh = Qh then gives 

Ul-QH)vof > \\lH{I-QH)Vof+Ul-lH)iI-QH)Vof 
= \\{lH-QH)Vof + \\{I-lH)Vof 

= wQnvor + wvor 



So 



Since 



and 



then 



-2 



v&Vn \\V\\a Folia ll^olla 



j=-n j=0 j=0 



^oWl = \\dxVof + 49yVof + e/h, 



K>h-^t€>h-^.h. ^ ' 



□ 



The above results show that in case of a grid that is aligned with the anisotropy direction 
the convergence of a standard two-level method (point-wise smoother and standard coars- 
ening) will deteriorate. One easily sees that for e < /i^ we get a poor convergence rate (no 
better than 1 - 0{h^)). 

However, the next result shows that when the grid is not aligned with the anisotropy 
direction (e.g., angle of rotation u = vr/4) the lower bound given in Lemma 3.1 does not 
apply and the standard two-level method is uniformly convergent in this case. 
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Lemma 3.3. Assume that is obtained from fio by a rotation with angle of rotation u = j. 
Then the error propagation operator corresponding to the two-level iteration with coarse space 
Vh and point-wise Gauss-Seidel smoother is a uniform contraction in the energy norm. In 
fact, we have the estimate 

(3.5) K = sup < C, 

vevi, \\v\\i 

and so 



(3.6) II^tlI 
with constant C independent of t and h. 



< 1 



1 



Proof. The proof follows the same lines as the proof of Theorem |3.4[ which will be given in 
the Section |6l □ 

Let us remark here that for decreasing values of the angle of rotation (i.e. decreasing uj 
from 7r/4 to 0) the convergence rate HE'r^.lU of a two-level method with point- wise smoother 
deteriorates as the angle of rotation becomes smaller. 

From the above considerations, it is clear that even in the case of aligned anisotropy 
one needs to use a special smoother or coarsening strategy in order to achieve uniform 
convergence. Our analysis shows that the line smoother or more generally a block smoother 
with blocks consisting of degrees of freedom along the anisotropy results in a uniformly 
convergent method. The Theorem below provides a uniform estimate on the convergence 
rate of the error propagation operator and is the main result in this paper. 

Theorem 3.4. For any angle of rotation uj G [0,7r], the two-level iteration with coarse space 
Vh and line (block) Gauss-Seidel smoother is a uniformly convergent method. In fact, we 
have 

(3.7) \\EtlV ' ^ 



< 1 



with constant C independent of e and h. 

The proof of this theorem is postponed to Section |6| The result follows from the stability 
and interpolation estimates that are given in Section |4] and Section [5] and Lemma 3^ 



4. Stability of the coarse grid interpolant 
In this section we prove the stability of the coarse grid interpolant. 

4 

In what follows, given a triangle K G Th {K = [J Ki with Ki G Th) as shown in the 

1=1 



Figure 4.1, we shall frequently use the following equalities 



2/2 - 


y& = 2/6 


- 2/1 


= 2/4 


- 2/5, 


2/3 - 


2/4 = 2/4 


- 2/2 


= 2/5 


- 2/6, 


yi - 


2/5 = 2/5 


- 2/1 


= 2/6 


- 2/4, 


\K\ 


= 4|ir/|. 


/ = 


= 1.. 


4. 



(4.1) 



Proposition 4.1. For any v G Vh, we have the following relation: 

3 



(4.2) 



dx 



K 




dv 
dx 
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Figure 4.1. A refinement of a coarse element K 



Proof. From (2.5) witli Ihv instead of f, we have 

d{lHv) 



dx 



K 



E&dK 



and from (2.5 ) with Ki instead of for / = 1, . . . , 4, we have 

dv 
dx 



K, 



EedKi 



Combine the above two equations and (4.1), it is immediate to verify the result. 



□ 



We are now ready to prove our first stability estimate. Since we have anisotropic diffusion 



problem in hand, we need to estimate separately 
the next Lemma. 

,, 2 

Lemma 4.2. For any v G Vh, we have 



dx 



and 



d{lHv) 



dx 



< 41 



dvu2 
dx WO' 



and 



dy 



, which is done in 



dy 



< 4 



I 9v 112 
IdyWo- 



Proof. We only need to prove this estimate locally for any K &Th- So we fix G Th and 



we would like to show that W^^W^k < 4||£llo,i^- 



From (4.2), for the L2 norm || ^^q^^^ IIo k have 



d{lHv) 



dx 



0,K 



< 




1=1 



Summing over all the elements then gives || IIo ^ ^llff IIo- -'■^ ^ similar fashion we can 

□ 



prove that ||^^^|^||o < 4|||^||o, and the proof of the lemma is complete. 

As a consequence, we have the following approximation result for the coarse grid inter- 
polant. 
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Lemma 4.3. For any v ^ Vh, we have 



d{v - Ihv) 


2 

< 


dv 


2 


d{v - Ihv) 


2 

< 


dv 


dx 





dx 





dy 





dy 



and 



lHv\\l<h'\v\l. 



Proof. The first two estimates follow from the inequalities given in Lemma 4^ The third 
estimate can be found in fi2[ Lemma 4.4]. □ 



Remark 4.4. In fact, in the proof of Lemma 4-^ there is no any requirement for the 
partition. The result is true for partition Th obtaining from regular refinement of any given 
partition Th- So is the Lemma \4.3[ 



5. Stability estimates on the fine grid 

In this section we give estimates on the stability of the partition of unity introduced in 
Section[2| equation (2.3). In what follows, to avoid proliferation of indicies, we will omit the 
subscript i and we will write 6 instead of 6i. 



4.1 



we label with 1, 2, 



For any given K G Th {K = [J Ki with Ki G Th) shown in Figure 

1=1 

and 3 the vertices of K, and 4, 5, and 6 the midpoints of K (these are also vertices of K4). 
The corresponding coordinates are denoted by {xj,yj), 1 < j < 6. Further, for a continuous 



function v, when there is no confusion, we write Vj :- 
In what follows we also denote 



'VK.Xj.yj, 



(5.1) 



Emin = arg min 



where 6^'^y is defined in (2.4). 



Then -Emm is related to the anisotropic direction. In fact, to indicate the dependence on the 
particular element, one may write -E^in instead of Emin, but for simplicity we have chosen 

K' 

to omit the superscript K4. Furthermore, we may denote -E^in = '^^5' {I^e*!/!}- 

Let us now consider a function w & Vh vanishing at the coarse grid vertices, that is, w 
satisfies Ihw = 0. We have the following 4 cases on a fixed K E Th- 

Case 0. 9 is zero in K4; 

Case 1. 6 is nonzero in K4 and convex in K4 (i.e. 6* 7^ a.e. in K4)] 
Case 2. 6 is nonzero at only one of the vertices of K4 and concave in K4; 
Case 3. 6 is nonzero at exactly two of the vertices of K4 and concave in K4. 

The rest of this section contains technical results and their proofs, which can be classified 
according to the cases above. To prove the stability estimates on the fine grid, we need to 
bound ll^^^^ll. 

• For the Case 0, there is nothing to prove, since in this case Ih{Ow) = 0. 



For the Case 1 the corresponding estimate is given in Lemma 5.2 In Case 1 we also 
need to assume quasi-uniformity of the mesh. We also note that Proposition 5.3 (for 
Case 1) contains an estimate which is later used in Case 2 and Case 3. 



The stability estimates in Case 2 and Case 3 are given in Lemmas |5.4| and |5.5 
respectively, under the assumption that Th is a. uniform partition. 
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Remark 5.1. In summary, for uniform mesh we have proved the stability estimate in all 
cases. In addition, we have proved some of the results in more general case of unstructured, 
but quasi-uniform mesh ( Case 1 ). 

Lemma 5.2. Assume that 6 ^ in K4 and convex in K4 (Case 1). Then 



d{h{ew)) 


2 

< 


d{9w) 


2 


dx 


0,K 


dx 


0,K 



Proof. Since Ihw = we obviously have that Inidw) = as well. From (4.2) in Proposition 
4.1 with V = IhiOw) we obtain that 



(5.2) 



E 

1=1 



d{h{ew)) 



dx 



d{h{ew)) 



K, 



dx 



d{iH{ih{ew))) 



K4 



dx 



d{iH{ew)) 



K 



dx 



K 



In addition, from (4.2), with v = w we have 
(5,3) 



E 

1=1 



dw 



dx 



K, 



dw 
dx 



MIhw) 



K4 



dx 



0. 



K 



Therefore, from (5.2) for the norm || ^^^A|^^||o,i<: we have: 



d{Ih{Ow)) 



dx 



0,K 



(5.4) 



1=1 
1 



dx 



Ei*-^ 



1=1 



d{h{dw)) 



dx 



{[(2/4 - VeWeWe + (1/5 - y4)05W5] 

+[{y5 - y4)0iWi + (?/6 - yb)0&WQf 
+[(z/6 - yb)d^w^ + (1/4 - y&iOiWif}. 



On the other hand, since ^ is a convex function in i^4, and 9 is supported in a 2h width 
strip, 9 should be convex in at least three of = 1 : 4). 
For any Ki in which 9 is convex, we have 



\e\\lK,>mY.9\xf,yf% 



where {xf'^,yf') denote the coordinates of the j-th vertex of element Ki for j = 1 : 3. Since 
the mesh is quasi-uniform, we have that max{yf'-} — min{?/f'} > h, and there exists at 

j=l:3 ■' i=l:3 ■' 

least one vertex jo such that 9{xjf^,yjQ) = 9{yjQ) > 1. Hence, if 6* is a convex function in 
Ki, then 116*110^^ > \Ki\, and we conclude that there are at least three elements Ki, where 



\9\\l^^ > \K,\ holds. 
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From this argument and (5.3) we get 



d{9w) 


2 


"1^ 


2 


dx 


0,K 


OX 


0,K 



K, 



1=1 



K, 



1=1 



K, 



(5.5) 



\K 



■{[(l/4 - y6)u!6 + (?/5 - ^4)^5]^ 

+[(2/5 - y4)w4 + ive - y5)we? 
+[(?/6 - y5)w5 + (1/4 - y6)w4]'^}. 



Introducing now 

M 



y5-y4 2/4 - 2/6 ' 

1/5-2/4 ye- ?/5 

,2/4 - 2/6 2/6 - 2/5 



e 



we rewrite (5.4) as 



d{h{9w)) 



dx 



0,K 



\K\ 



\MQz\ 



while (5.5) can be rewritten as 

d{ew 



dx 



0,K 



\K\ 



\Mz\ 



^2' 



{w4,W5,Wq) 



{W4,W5,W6) 



Here, || • We^ is the usual Euclidean norm on M^. 

To prove the estimate || 9{ih{dw)) ||2 ^ ^ || ^^"'^ ||o o^^Y need to show that 



(5.6) 



dx 

^ "MQzWl < ^ 



\K 



K\ 



\Mz\ 



for all z G 



Such an inequality is easy to get in the case of det(M) = 0, so we may assume M 
is invertible (i.e. det(M) = 2(^/6 — 2/5) (2/4 ~ 2/6) (2/5 ~ 2/4) 0)- We then need a bound 
on the eigenvalues of M'^QM'^eM-^ = {M-^QM){M-^Q)MY . In fact, we only need to 
bound the entries of M~^QM because all the norms of this 3x3 matrix are equivalent. 
Thus, if the entries of M~^QM are bounded in absolute value, then the eigenvalues of 
{M~^QM){M~^QMy are bounded and consequently (5.6) holds. 

Directly computing the inverse of M gives 



det(M) 



-(2/6-2/5)^ 
(2/4 - 2/6) (2/6 - 2/5) 



(2/6 - 2/5) (2/4 - 2/6) {y& - yb){y^ - yi) 



■(2/4-2/6) 



(2/5 - 2/4) (2/6 - 2/5) (2/5 - 2/4) (2/4 - 2/6) 



(2/4 - 2/6) (2/5 - 2/4) 
-(2/5-2/4)^ 



We then calculate M ^BM to obtain that 



^5 + 



M-^QM 



VA-yg I 



— ^'5 + ^'g 



— 6'4 + 6'6 



^4 + Or 



-^6 



S/6-S/5 ^ S/4-J/6 

Since 9 is convex in i^^4, by the definition of 9, it is easy to see that 

l^'e - 9^\ < h'^lye - yg], 
\95 - 94\ < h'^lyr, - ?/4|, 



—94 + 9^ 
(94 - 95) 
94 + 9^ 



ys-yi ^ 
y4-y6 
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Since \yi — yj\ ^ we have \{M~^QM)ij\ < 1 and the proof of the Lemma is complete. □ 



Next result is an auxiliary estimate used later in the proof of Lemma |5.4| and 5^ In the 
statement of the lemma we used the notation given at the end of Section |2} 

Proposition 5.3. Assume that 9^0 and convex in K4^. Then the following inequality holds 



d{9w) 



dx 



0,K 



> \K\{max{\6^^y\r ■ KlJ^ + ■ max 



Proof. Let Emin be defined as (5.1 ). Without loss of generality, assume Emin = {(a^4, Va), {x^, 1/5)}, 
and then wq = Wp^ . This means 



bs - I/4I = min{|?/5 - y^l, |?/4 - y^l, \ye - y^l}. 



Hence 



Vb - 


2/4 


<1, 


2/5 


-2/4 


<1, 


ye - 


2/5 




2/4 


-2/6 


have 












2/4 - 


2/6 


<2, 


2/6 


-2/5 


< 2. 


2/6 - 


2/5 




2/4 


-2/6 





According to the expression (5.5) and the above inequalities, we have 
2 1 



d{9w) 



dx 



0,K 



\K 



{[(2/4 - 2/6)^6 + (2/5 - 2/4)^5]^ + [(2/5 - 2/4)^4 + (2/6 - 2/5)^^6]' 



> 



+ [(2/6 - 2/5)^5 + (2/4 - 2/6)^^4]^} 

{[(2/4 - y&)'WQ + (^5 - 2/4)^5] + [(?/5 - y4)ti'4 + (2/6 - 2/5)^^6] 



yA-y& 
2/6 - 2/5 



-[(2/6 - 2/5)^5 + (2/4 - 2/6)^^4] 



2/5-2/4 -12 
2/6 - 2/5 



> 



2 

1^ 
1 



[(2/4 - 2/6)w^6]^ 

max{|?/5 - ?/4|, |?/4 - yel, I2/6 - yb\Ywl. 



Combining with (5.5), we have 

2 



d{9w) 



dx 



1 



> r^{[(2/5 - 2/4)^4]^ + [(2/5 - 2/4)^5]^}. 



So 

d{9w) 



dx 



0,K 



~ jKl^^^^^^y^ ~ ^^1' ~ ^el, I2/6 - 2/51)^^6 + (1/5 - 2/4)^^4 + (2/5 - yifwl}. 



Notice again, here Emin = {(x4,?/4), (xs, 2/5)} and we = . , then we get the result. □ 
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We need to notice that till now we only require the mesh is quasi-uniform, since when 6 is 
convex in element K4, the semi-norm of interpolation function ||^^^^^|^^||o,x can be bounded 

by II ^gj^ ||o,i<:- However, this is not true when 9{y) is concave. In this case, II ^''^'gf"''*'' llo,J<: 
may also depend on some neighboring element. To get the information of the neighboring 
element, we assume the partition 7h is uniform in the following. 

Lemma 5.4. Assume that 9 is nonzero at only one vertex of K4 and that K' is the unique 
element from Th which has this vertex on one of its edges. Assume also that 9 is concave in 
Ki (Case 2). Then the following inequality holds 



d{Ih{9w)) 


2 

< 


d{9w) 


2 


dx 


0,K 


dx 


0,K 




Figure 5.1. The coarse elements K and K' sharing the same midpoint 



Proof. Without loss of generality, assume 9^^ is the only nonzero value. There are two 
possibilities: (a) Eq = Emin] and (b) Eq ^ Emin- 

Proof in case (a). Since Eq = Emm, we conclude that |5^gl/| 
We then have 

\K\ 



min {\6^^y\}. 



d{h{0w)) 



dx 



0,K 



< 



\K\ max {\6^^y\} 



Since the partition Th is uniform (see Figure 5.1), and K' is the element sharing the same 
point (a^Eo'^Bo) with K, we know that the values of 9 at midpoints of K' are all nonzero. 
This is so, because the support of 9, whose width is 2h must include K' in its interior. Assume 



now that Eq is the edge opposite to point {x^^,y^^) in K'^ (i.e. {x^^^,y^^^ 
Figure 



5.1). Observe that E'q = E'^ 



I 



min {l^g/'yl}, because E'q is a parallel 



^4 K\ 
^e'.^Ve'J^ 



see 



E'edK: 
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translation of Eq. By Proposition 5^, we now have 

2 



d{ew) 



dx 



0,K' 



> \K'\ maxj\6^h\nw^4r = 1^1 max 



So II ^^^'dx^^^ llo.A" ~ II ^aT^ llo E"'' and this completes the proof in case (a). 
Proof in case (b). In case (b) we have Eq ^ Emm and hence |5^q*2/| 7^ niin {|5^''y|}- 

'Eq 



E&dKi 



Since 9^f^ is the only nonzero value among the values of 9 at the vertices of K/^, we easily get 



Then 



d{h{ew)) 



dx 



0,K 



\K\ 



< 



E \^E'y\'i(^E:^E:r (irom^) 

EGdK4 



1^1 max{|<5^^y|}^(^^^^)^«^ 



< \K\ max {2\K,\5fy\Yh-' min {|5^^1/|}^«: 



E&dKi 

< \K\mm{\6j^^y\nw^^r. 



E&dKi 



Let K' be the same as before, then by Proposition 5.3, we have 

2 



d{ew) 



dx 



0,K' 



> |ir'|^min,{|4^y|}2(i.Jf)2 = minj|5^^y|}2«^ 



E€dK4 



Combining the last two inequalities then gives || ^^^'^q^^^^ WI^k ~ II '^aT'' llo x'- This completes 
the proof in case (b), and also the proof of the Lemma. □ 

The next Lemma gives the stability estimates in the last case (Case 3) and we refer to 
Figure 5.2 for clarifying the notation. 

Lemma 5.5. Assume that 9 is nonzero at exactly two vertices of and concave in 
(Case 3). Let K' be an element from Th which shares with K^^ the vertex at which 9 has 
larger value on one of its edges. Then the following inequality holds 



d{h{9w)) 


2 

< 


d{9w) 


2 


dx 


0,K 


dx 


0,K[JK 



Proof. Without loss of generality we may assume that ^4 and 9^ are the only nonzero values 
of 9, we may also assume that 9^ > 9^. As a consequence, K' will share (x4,?/4) with K. 
We consider two possibilities: (a) Emin = {(2:4,1/4), (3:5, Z/5)}; (b) -E^m 7^ {(x4,i/4), (a;5,l/5)}- 
Proof of (a). Since E^in = {(0:4, 1/4), (3:5,2/5)}, we have that I2/5-2/4I = ram{\y5-y^\, ly^- 
Vel, Ive ~ Ubl}- Hence, from (2.5) we obtain 



dilhi9w)) 



dx 



0,K 



\K 



■{((2/5 - yA)9iWif + ((?/5 - yi)9^w^f 
+((2/6 - yb)9^w^ + (2/4 - y&)9AWif}. 
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Figure 5.2. The coarse elements K and K' sharing the same midpoint 



We want to bound now all the terms on the right side of the above relation with quantities 
independent of the values of 9. From the fact that 6* < 1 we have that 

The other terms are bounded as follows 



iiVb - y4:)05W5f < ((^5 - Vijw^f < ((ys - y4)w5) 



2/1/6 -1/5x2 



yi-y& 



< 



{{{ye - yb)wb + (2/4 - y&)wA)— — — - {y^ - yA)wAf 

yA - ye 

((Z/6 - y5)w5 + (?/4 - y&)Wif{— —f + ((?/5 - yA)wAf 

yA-y& 



^ ((?/6 - yb)w^ + (1/4 - y&)wi) + ((z/5 - yA)wi) , 

and also 

((2/6 - ybW^w^ + (2/4 - 2/6)6*4^4)^ = ((2/6 - ybW^w^ + (1/4 - y6)05W4 - (1/4 - ye) (6*5 - ^'4)^4)^ 

= (((2/6 - 2/5)^5 + (2/4 - 2/6)w4)6'5 - (2/4 - y&)wJ^^-T^f 



h 



Hence we get 



dx 



0,K 



^ ((ye - y5)w5 + (y4 - y6)w4f + ((2/5 - ^4)^4)^. 



~ r^^((^5 - 2/4)^^4) + ((2/6 - 2/5)^^5 + (2/4 - 2/6)^4) }. 



In this case, 



d{ew) 



dx 



0,K 



~ r^^((^6 - 2/5)^5 + (2/4 - 2/6)^^4) }■ 
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Since K' denotes the element sharing (x4,?/4) with K, we obtain 
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d{ew) 



dx 



1 



.K'\ i„.K' 



K'l i„.K' 



.K'l^2„,,2^ 



0,K' 



{min{|y5 - y^l, {y^ - ye\, \ye - y^]} w^} 



so 



d{ew) 



dx 



0,K[JK' 



1 

T^Aiy^-yifwl}, 



~ r^^((^5 - 2/4)^4)^ + {{ye - y5)w5 + (2/4 - 1/6)^4)^} > 



dx 



0,K 



This completes the proof in case (a). 

Proof of (b). In this case we have that Emm 7^ {(x4,?/4), (x5,?/5)}, which is equivalent to 
I1/5 - 2/4I > min{|y5 - y^l, {y^ - ye\, {y^ - y^]}. 

Since 6*4 > 9^ (6*4 = 9^ can not be true in this case, because 6*4 = 9^ implies y^ = y^), we 
can get \ye-y5\ < lyA-yel- It is then easy to see \ye-y5\ = niin{|?/5-?/4|, |?/4-?/6|, |?/6-?/5|}. 



Then 



9^<h ^ min{|y5 - t/4|, |?/4 - y&\Ay& - y-^W = h ^\yQ - ys]. 



So from (2.5), we have 
2 



d{h{9w)) 



dx 



0,K 



< 



< 



1 

Jk 


{((2/5 


- y4)95W5f + ((?/5 


- ?/4)6'4ti'4)^ + {{ye 


- 2/5)6*5^5 + (2/4 - y6)94W4f} 


1 

W 


{((2/5 


- y4)d5W5)'^ + ((?/5 


- y4)94W4f + {{ye 


- yb)d5W^f + ((2/4 - 2/6)6'4W4)^} 


1 

w 


{((1/6 


- Z/5)W^5)^ + ((1/4 - 


y6)w4f}. 





In this case, 



d{9w) 



dx 



0,K 



~ r^^((^6 - ^5)^5 + (?/4 - 1/6)^4) }• 



Since the partition Th is uniform (see Figure 5.2), and K' is the element sharing the same 
point (x4,?/4) with K, we know that the values of 9 at midpoints of K' are all nonzero. 
Observe that the edge opposite to point (X4, y^) in K'^ is a parallel translation of the edge 
opposite to point (X4, y^) in K4. That is to say, point (X4, y^) in i^' (also point 4 in /T) is just 
the midpoint opposite to edge with mm{\y^ — 2/4" |, Iz/f —yQlAVe "Z/s' I}- Proposition 



5.3, we have 



d{9w) 



dx 



> ^{max{|^f - yf'l {yf - y-' {, |yf - 1}^} 



.K'l L.X 



.K'l^2„,,2^ 



0,K' 



\K 
1 

w 



{max{|?/5 - ?/4|, |?/4 - T/el, be - Z/sl} ^^4} 



Combining the last two inequalities, we have 
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d{ew) 



dx 



0,K\JK' 



> (^{((z/e - y5)w5 + {Vi - y6)wiY + max{|?/5 - y^l, {y^ - y^l, {y^ - y5\Yw 



> 



d{h{ew)) 



dx 



This completes the proof of case (b), and also the proof of Lemma. 



Lemma 5.6. For any w G Vh, if IhW = 0, then for any 1 < i < L, 



□ 



dihieiw)) 



dx 



< 



dx 



Proof. The estimate follows from the local (element-wise) estimates given by Lemma 5.2 
5.4[ 5.5 and summation over all elements from Th- D 



6. Proof of the theorem 13.41 
In this section we prove the convergence result that we have already stated in Section [3j 



Theorem |3.4[ For any angle of rotation lo G [0,7r], the two-level iteration with coarse space 
Vh and line (block) Gauss-Seidel smoother is a uniformly convergent method. In fact, we 
have 

1 



\E' 



TL\\a 



< 1 



with constant C independent of e and h. 
Proof. From Lemma 



3.2 



we only need to prove sup ^^'^^ imi2^"^^* ^ C. 



For any v G Vh, let w := v — Ihv, Wi := lh{0i{y)w), it is easy to see Ihw = 0, and 
^Wi = ^Ih{Oi{y)w) = Ih^iOi{y)w) = hw 



w. 



Then 



vev^ Ml v&v^ Wnl 



\w\ 



sup " = sup inf 



< 



sup 



So we only need to prove for any v G Vh, MiWa 
First, the decomposition is stable in L^, 



< \\vP. 
~ II \\a 



(6.1) 



\w 



I? < 
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Since the sum is along x direction, 

i i 

^ ^\\9x{0i{y)w)\\l (by Lemma [K6|) 

i 
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(6.2 



then 



< 



< 



2|| o ||2 



E^ ll^^'^illo + E/ ^^"^11^*11 inverse inequahty) 



eh-^\\w\\f, (by (6.1) and (6.2)) 



\d^{v - Ihv)\\(^ + eh \\v - Ihv 



^ \\d^v\\Q + eh'^h^\v\i (by Lemma 4.3) 



□ 



7. Numerical Experiments 

7.1. Tests for two-level method on a rotated uniform mesh. We first test the perfor- 
mance of the two-level iterative method and its convergence properties with respect to e and 
h. We pick as initial triangulation a 4 x 4 mesh with a characteristic mesh size ho = \\/2 
as shown in Figure 2J^ We then apply the two-level method described earlier on sequence 
of meshes with mesh sizes h^ = 2~^hQ, = 1, . . . , 6. 



The energy norm of the error of two- level method is depicted in Figure 7.1 show 



that two-level method is uniformly convergent w.r.t. e and h, which agree with the theoretical 
results we have proved in the previous sections. 



7.2. Tests for two-level method on a general unstructured mesh. Similarly to the 
case of uniform mesh, for a general unstructured mesh we choose h^ = 0.9 as the maximum 
diameter of the triangles on the coarsest mesh To as shown in Figure 7.2 This coarsest 



mesh is then refined 6 times and get the mesh to obtain a sequence of triangulations with 
characteristic mesh sizes h = 2~''hQ, k = 1, . . . ,6. 

The energy norm of the error propagation operator for the two level method, are 
shown in Figure 7.3 The uniform convergence is clearly seen from the plots. Theoretical 



justification of such uniform convergence is however much more difficult and is a topic of 
current and future research. 
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Figure 7.2. Plot of unstructured grids used in the numerical examples for three 
values of the angle of rotation of anisotropy (X' = 0, (X' = 7r/6 and uo = tt/A. 
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